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a natural family of compatible Poisson structures for the Hamiltonian formulation of 
Lax equations. Examples for which our formalism applies include the Benny hier- 
. achy, the dispersionless Toda lattice hierachy, the dispersionless KP and modified KP 

f*^*) hierachies, the dispersionless Dym hierachy etc. 
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i 1. Introduction. 

Two Poisson brackets on the same manifold are said to be compatible if their 
sum is also a Poisson bracket [GDO, M]. There are many examples of integrable 
^ ■ systems which are Hamiltonian with respect to two compatible Poisson structures 
(see, e.g. [DO]). Indeed, when one of the structures happens to be nondegenerate, 
there is a simple way which allows one to produce a whole family of compatible 
Poisson structures [KR, RSTS1]. However, the existence of further structures is not 
a necessity when the two compatible structures are both degenerate. 

In the late seventies, we saw the beginning of the Lie algebraic approach to in- 
tegrable systems [K, A]. The Korteweg de-Vries (KdV) equation, for example, was 
shown to be a Hamiltonian system on coadjoint orbits [A] . Furthermore, the second 
Poisson structure for KdV type equations was constructed on subspaces of the al- 
gebra of formal pseudo-differential operators [A, GD]. We now refer to this second 
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structure as the Adler-Gelfand-Dickey structure. Recently, it was found to be of 
independent interest in conformal field theory [DFIZ]. In the mean time, the Lie 
algebraic approach to integrable systems was extensively developed, particularly by 
the Russian school in St. Petersburg (see, e.g., the survey in [RSTS2]). In the 
so-called r-matrix framework, the simplest Poisson structures for the Hamilton- 
ian formulation of Lax equations on Lie algebras are the linear Poisson structures 
associated with the /^-brackets. In the case where q is the Lie algebra of a non- 
commutative, associative algebra, a construction of quadratic brackets which give 
Lax equations was first available for the skew-symmetric r-matrices satisfying the 
modified Yang-Baxter equation [STS1]. Subsequently, this was superseded by a 
more general construction valid for a wider class of r-matrices [LP1,LP2]. Indeed, 
in [LP2] , even a third order structure was found. At this juncture, the reader should 
note that on the abstract level of associative algebras, neither the linear structure 
nor the quadratic structure is nondegenerate. Therefore, the recipe for producing a 
whole family of structures is not applicable in this context. As a matter of fact, no 
Poisson structures with order > 3 was ever found. In this connection, we would like 
to mention the thesis of Strack [ST] , which showed (by using computer algebra) that 
beyond order 3, no Poisson structures of a certain form can exist for the Hamilton- 
ian formulation of Lax equations. So this is the state of affairs for noncommutative, 
associative algebras. 

In this paper, we address the Hamiltonian formulation of Lax equations, as be- 
fore, but in the context of Poisson algebras. Here, we show how to construct a 
natural family of compatible Poisson structures on the full algebra. On the group 
of invertible elements (if non-empty and forms an open subset), similar considera- 
tion shows we can even define structures of negative order. Thus the situation for 
Poisson algebras, in which multiplication is commutative, is entirely different. Re- 
call that a Poisson algebra is by definition a commutative, associative algebra with 
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unit 1 equipped with a Lie bracket such that the Leibniz rule holds [Wl] . The most 
familiar examples of Poisson algebras are given by the collection of smooth functions 
on Poisson manifolds. For us, the particular examples which have partly motivated 
this work are the algebras associated with the truncated Benney's equation [G- 
KR], and the various dispersionless equations [DM, K, TT] which are currently of 
interest in topological field theory [D, K]. As the reader will see, a family of vector 
fields V n , n > — 1, plays the key role in this investigation. These vector fields V n 
are invariants of degree 1 of the vector fields associated with the Lax equations, 
and satisfy the Virasoro relations [V m , V n ] = (n — m)V m+n . For a given classical 
r- matrix on the Poisson algebra, we can construct the associated linear bracket. 
If we denote by 7r_i the bivector field corresponding to this basic linear structure, 
we shall show that the Lie derivatives L Vm K-i essentially generate all higher order 
structures. Thus our construction works for an arbitrary classical r-matrix! This 
is in marked contrast to previous results on quadratic Poisson structures on non- 
commutative, associative algebras [STS1, LP1, LP2], where one has to make rather 
stringent assumptions on the r-matrix. In this connection, we would like to remind 
the reader of the important difference between the notions of double Lie algebras 
and Lie bialgebras. Recall that the former was motivated by the study of integrable 
systems [STS1] and is associated with classical r-matrices. On the other hand, the 
notion of Lie bialgebras had its origin in the geometry of Poisson Lie groups [DR] . 
The two do intersect, for example, in the class of double Lie algebras called Baxter 
Lie algebras [STS2] (where the r-matrix satisfies additional properties) . In our case, 
as the r-matrix is assumed to be completely arbitrary, we are working within the 
framework of double Lie algebras here. 

The paper is organized as follows. In sec. 2, we assemble a number of basic facts 
and definitions which will be used in the paper. In Sec. 3, we formulate the main 
result and display the explicit formulas for the linear, quadratic, and higher order 



4 LUEN-CHAU LI 

structures. Then we study a number of basic properties. In order to prepare for the 
proof of the main result, we introduce the vector fields V n in Sect. 4 and discuss 
their relation with the Lax equations. Then, in Sect. 5, we give a proof of the main 
result. In order to illustrate the use of our construction in Sect. 3, we describe the 
multi-Hamiltonian formalism of some concrete partial differential equations in Sect. 
6. Our examples include the hierachy of truncated Benny equations [B, G-KR] in 
nonlinear waves, the dispersionless Toda lattice hierachy [DM], the dispersionless 
KP [K, TT] and modified KP hierachies, and the dispersionless Dym hierachy. Note 
that in each example, the set of Lax operators under consideration is a submanifold 
of the full Poisson algebra. However, this submanifold is not necessarily a Poisson 
submanifold of the full algebra equipped with a bracket which comes from Sect. 3. 
For this reason, the passage from the bracket on the algebra to the Hamiltonian 
structure on the submanifold of Lax operators might involve the process of reduction 
[MR]. Thus in our examples, we find Dirac reduction [D, MR] (i.e. reduction with 
constraints) comes in naturally. For the Benny hierachy and the dispersionless Toda 
lattice hierachy, we shall compute the first few Poisson structures explicitly, and 
illustrate the use of Dirac reduction. Our explicit expressions for the structures not 
only allow us to find the Casimir functions, they also show that the structures which 
come from our Poisson algebras are of hydrodynamic type or its generalizations [DN, 
F]. Indeed, as it turns out, all the higher structures of the dispersionless Toda lattice 
hierachy are nonlocal generalizations of brackets of hydrodynamic type. This shows 
how our construction in Sect. 3 can get complicated upon reduction to a specific 
submanifold of Lax operators. 

To close we stress again that our main result is formulated along the lines of 
the r-matrix approach (where Poisson structures are defined either on Lie algebras 
or their duals, or on Lie groups) and applies to all Poisson algebras satisfying the 
assumptions of Theorem 3.2. In any concrete applications, the use of reduction 
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techniques (where necessary) is perfectly natural and the reader should not feel 
uncomfortable under such circumstances. 

2. Preliminaries. 

We collect in this section a number of basic facts, and introduce some terminology 
which will be used in the sequel. 

Let P be a smooth manifold. A Poisson bracket {•, •} on P is a Lie bracket on 
C°°(P) which satisfies the derivation property in each argument. If it is the bivector 
field corresponding to the bracket operation, i.e. 

(2.1) {F, H} = n(dF, dH) , 

then it is well-known that the Jacobi identity for {■, •} is equivalent to [7r,7r]s = 
[W2], where [•, -}s is the Schouten bracket [S]. Recall that if T(A k TM) is the space 
of sections of the vector bundle A fc TM, and A*(M) = © T(A fc TM), the Schouten 

fc>0 

bracket [■, -]s is the bilinear map 

(2.2) [;-] s : A*(M) x A*(M) -> A*(M) 

which extends the usual Lie bracket operation on T(TM) and makes A*(M) into a 
Lie superalgebra. In particular, the following graded Jacobi identity holds: 

(2.3) (-l) p >, [v, w]s]s + (-l) qp [v, [w, u] s ]s + (-l) r >, [u, v] s ]s = 

where u E T(A P TM), v E T(A q TM) and w E T(A r TM). 

As we mentioned in the introduction, two Poisson brackets on P are said to be 
compatible if their sum is also a Poisson bracket, i.e. satisfies the Jacobi identity 
[GDO, M]. In terms of the corresponding bivector fields 7Ti and 7r 2 , this is equivalent 
to [tti, tv 2 ]s = 0, as [7Tj, TVi]s = 0, % = 1, 2. 

In this paper, we shall construct compatible Poisson structures for the Hamilton- 
ian formulation of Lax equations (associated with r-matrices) when the underlying 
manifold P is a Poisson algebra. 



6 



LUEN-CHAU LI 



Definition 2.4. Let A be a commutative, associative algebra with unit 1. If there 
is a Lie bracket on A such that for each element a G A, the operator ad a : b i— > [a, 6] 
is a derivation of the multiplication, then (A, [•,•]) is called a Poisson algebra. 

Thus the Poisson algebras are Lie algebras with an additional associative algebra 
structure (with commutative multiplication and unit 1) related by the derivation 
property to the Lie bracket. Note that some authors call the Lie bracket on A the 
Poisson structure on A (see, for example, [Wl]), but we shall refrain from such 
usage in order to avoid confusion. 

We now recall the notion of a classical r-matrix [STS1]. Let q be a Lie algebra. 
A linear operator R in the space q is called a classical r-matrix if the i?-bracket 
given by 



is a Lie bracket, i.e. satisifes the Jacobi identity. Some well-known sufficient condi- 
tions for R G End(g) to be a classical r-matrix are the Yang-Baxter equation and 
the modified Yang-Baxter equation. But in this paper, we can establish our results 
without assuming these conditions. 

To close this section, we define what we mean by Lax equations. 

Definition 2.6. Let A be a Poisson algebra, and suppose R G End(A) is a classical 
r-matrix. Equations of the form 



(2.5) 



[X,Y] R 



±([RX,Y] + [X,RY]), 



X 7 Yeg 



(2.7) 



L=[R(X(L)),L], 



LeA 



where X : A — > A is a smooth map satisfying 



(2.8) 



[X(L),L} = 0, 



dX(L).[L',L] = [L',X(L)}, 



L,L' <E A 



are called Lax equations. 
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The basic Lax equations on A are given by 

(2.9) L = Z m (L) = [R(L m ),Ll m>l. 

More generally, if H is a smooth ad-invariant function (in the sense defined in (3.1)), 
then L — [R(L m dH(L)), L] is also a Lax equation. 

3. A Family of Compatible Poisson Structures on Poisson Algebras. 

In what follows, we shall assume the Poisson algebra A is equipped with a non- 
degenerate ad-invariant pairing (•, •). A function F defined on A is said to be smooth 
if there exists a map dF : A — > A such that 

(3.1) ±\ t=0 F(L + tL') = (dF(L),L') , L,L'eA 

Theorem 3.2. Let A be a Poisson algebra with Lie bracket [•, •] and non- degenerate 
ad-invariant pairing (•, •) with respect to which the operation of multiplication is 
symmetric, i.e. (XY, Z) = (X,YZ), VX,Y,Z e A. Assume R G End(A) is a 
classical r-matrix, then 

(a) for each integer n > — 1, the formula 

(3.3) {F,H} {n) (L) = (L, [R(L n+1 dF(L)),dH(L)} + [dF(L), R(L n+1 dH(L))]) 

(where F and H are smooth) defines a Poisson structure on A, 

(b) the structures {-,-}( n ) are compatible with each other, 

(c) if'Kn is the bivector field corresponding to {•, and D 1Tn : A* (A) — > A* (A) 
is the associated coboundary operator, i.e. D 7Vn X = [ir n ,X]s, X e A* (A). 
There exists vector fields V m on A, m > —1 satisfying the Virasoro relations 
[V m , V n ] = (n - m)V m+n such that D^Vm = (n - m)7i m+n , m,n> -1. 

□ 

We shall prove this result in Section 5, after we introduce the vector fields V m in 
Section 4 and explain what they are in relation to the Lax equations. As the reader 
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will see, the relations [7r n , V m ]s = (n — m)iT m+n between the bivector fields which 
we establish at the beginning of Section 5 play the key role in proving parts (a) and 
(b) of the above theorem. They are also responsible for the following. 

Corollary 3.4. (Involution of Casimir Functions) {H® n (A), H® n (A)}( m+n j = 0, 
m, n > — 1. m ^ n. 

Proof. This follows from the formula [7r n , V m ]s(dF, dH) = Lv m TT n (dF,dH) = 
V m {F, H} {n) - {V m F, H} (n) - {F, V m H} (n) . □ 

Remark 3. 5. Note that from the compatibility of the structures, it follows that 

(3-6) {H°JA),H°JA)} (n) <zH°JA). 

We now give a number of basic properties of the Poisson structures {•,•}(«)> 
n > -1. 

Theorem 3.7. (a) Smooth functions in A which are ad-invariant Poisson com- 
mute in {•, -}( n )- 

(b) The Hamiltonian system generated by a smooth ad-invariant function H in 
the Poisson structure {■, -}( n ) is given by the Lax equation L = [R(L n+1 dH(L)), L\. 

Proof, (a) If F and H are smooth functions in A which are ad-invariant, we have 
[dF(L),L] = [dH(L),L] = 0. Therefore, {F,H} {n) (L) = ([dH(L), L], R(L n+1 dF(L))) + 
([L,dF(L)}, R(L n+1 dH(L))) = 0. 

(b) If H is ad-invariant, for any smooth F, we have {F, H}^(L) = (L, [dF(L) : 
R(L n+1 dH(L))}) = (dF(L), [R(L n+1 dH(L)),L]). □ 

^From formula (3.3), it is clear that the bracket {•, -}( n ) vanishes at the unit 
1. Therefore, the linearization of {•,•}(„) defines a Lie bracket on A, and an easy 
calculation shows it coincides with the .R-bracket [•, •]#. 

The following result is reminiscent of the multiplicative property of Poisson Lie 
groups [DR]. However, it is in the context of a Poisson algebra and the reason for 
its validity is entirely different. 
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Theorem 3.8. Equip A with the structure {•, -}(o) o,nd A x A with the product 
structure. Then the multiplication map m : A x A — > A is a Poisson map. 

Proof. Let F and H be smooth functions on A. For Li, L 2 G A, let L = m(Li, L 2 ). 
Clearly, F o m depends on two variables and by taking its derivative with respect 
to the i-th variable, i = 1,2, we obtain d\{F o m)(L l5 L 2 ) = L 2 dF(L), d 2 (F o 
m)(L u L 2 ) = LidF(L). To simplify notation, let X x = dF(L), X 2 = dH(L) and 
denote the product structure on A x A also by {•, -}(o), then we have 

(*) {Fom,tf om} (0 )(Li,L 2 ) 

= (L u [#(LXi), L 2 X 2 ] + [L 2 X U R(LX 2 )]) 
+ (L 2 , [R(LXi), L1X2] + [L 1 X 1 , R(LX 2 )]) . 

By the derivation property of [•,•], the commutativity of multiplication and its 
symmetry with respect to the ad-invariant pairing (•, •), we have 
(L 1 ,[R(LX 1 ),L 2 X 2 ]) 

= (L, [i?(LX!),X 2 ]) - (L 2 , [R(LX 1 ),L 1 X 2 ]) . 

Likewise, 

(L 2 ,[L 1 X 1 ,R(LX 2 )]) 

= (L, [X u R(LX 2 )} - (Lx, [L 2 X U R(LX 2 )}) . 
When we insert these relations in (*), the result follows. □ 

Consider now A inv , the group of invertible elements of A. We assume A inv 7^ 
and form an open subset of A. Then we can define vector fields -Z_ m , V_ n for 
m > 1, n > 2, on A inv as in formulas (4.2) and (4.5). If we define 
(3.9) 

{F, H} { _ n) (L) = (L, [R(L~ n+1 dF(L)), dH(L)] + [dF(L), R(L~ n+1 dH(L))]), n > 2 

for smooth functions F and H on A inv , it is easy to check that the analysis in 
Section 5 also holds for these objects. In particular, this means {•, -}(- n ) are Poisson 
structures on A- inv . 
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Theorem 3.10. Let i : A inv — > A inv be the inversion map, i.e. t(L) = L _1 . Then 
{F o H o L\( n ) (L) = —{F, o i(L), n > 0, for all smooth functions F and H 

on A inv . 

Proof. We have d(F o t){L) = -L" 2 rfF(L" 1 ) and so 

{F o l,H o /,}( n )(L) = (L, [R(L n ~ 1 dF(L~ 1 )), L~ 2 dH "(L" 1 )] - (F <-> #)) . Now, 

(L, [i?(L"- 1 rfF(L- 1 )),L- 2 ^(L- 1 )]) 
= (L, L- 2 [i?(L n - 1 rfF(L- 1 )),rfif(L- 1 )]) + (L^(L- 1 ), [i?(L n_1 <iF(L~ 1 )), L" 2 ]) 
= (L" 1 , [i?(L n_1 (iF(L~ 1 )), dH(L- 1 )]) + 2(dH(L- 1 ), ^(L^F^- 1 )), L" 1 ]) 
= - (L-\ ^(^-^(L- 1 )), ^(L" 1 )]) . 

Hence the assertion follows. □ 

4. Lax Equations on Poisson Algebras and Virasoro Invariants. 

According to Definition 2.6, corresponding to each smooth map A : A — > A 
satisfying (2.8) is a Lax equation 

(4.1) L = A(L) = [ J R(A(L)),L] 

To prepare for the proof of Theorem 3.2, we shall introduce vector fields V n , n > — 1 
on A which are related to the Lax equations. Before we do so, we first prove 

Theorem 4.2. Let A, Y : A — > A 6e smooth maps satisfying (2.8). Then [A, F] = 
0. 



Proof. We have 



dF(L) • A(L) 

=[i2(dF(L) • A(L)), L] + [R(Y(L)), A(L)] 

= [fl([fl(X(L)), r(L)]),L] + [ J R(y(L)), [#(A(L)),L]] 
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Therefore, 
[X,Y](L) 

=2[R([X(L),Y(L)] R ), L] + [R(Y(L)), [R(X(L)),L]]-[R(X(L)), [R(Y(L)),L]] 
= [2R([X(L),Y(L)] R ), L] - [[R(X(L)), R(Y(L))], L] , by Jacobi identity 
= - [[R(X(L)), R(Y(L))} - 2R([X(L), Y(L)] R ), L] . 

Let B R (X,Y) = [RX, RY] - 2R([X,Y] R ). Then R is a classical r-matrix iff 
[B R (X,Y),Z] + [B R (Y,Z),X] + [B R (Z,X),Y] = 0, V X,Y,Z e A. Using the 
ad-invariant pairing, this is equivalent to 

[B R (X, Y), Z] =R*[RX, [Y, Z\\ - R*[X, R*[Y, Z\\ - [RX, R*[Y, Z\] 

+ R* [RY, [Z, X]]-R* [Y, R* [Z, X}} - [RY, R* [Z, X}} 

If we now put X = X(L), Y = Y(L) and Z = L in the above relation, we obtain 
[X, Y](L) = 0, as asserted. □ 

The vector fields V n , n > —1, are defined as follows: 
(4.3) V n (L) = L n+1 , n>-l 

Theorem 4.4. The vector fields V n satisfy the Virasoro relations [V m , V n ] = (n — 
m)V m+n , m,n> -1. 

Proof. Clear. □ 

Given a smooth manifold M and a vector field V on M, recall that a tensor field 
T is an invariant tensor field of V iff LyT = 0. Generalizing one step further, we 
shall say that T is an invariant tensor field of degree 1 iff L 2 X T = 0. The vector 
fields V n introduced in (4.3) above are invariants of degree 1 of the vector fields X 
corresponding to the Lax equations. Indeed, we have 

Theorem 4.5. If X : A — > A is a smooth map satisfying (2.8), we have Ly m X = Y 
where Y(L) = dX(L) ■ V m (L). 
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Proof. 

[V m ,X](L) 
=dX(L).V m (L)-dV m (L),X(L) 

= [R(dX(L) ■ V m (L)),L] + [R(X(L)),V m (L)} - (m + l)L m [R(X(L)), L] 

= [R(dX(L),V m (L)),L}. 
Thus, it remains to show Y(L) = dX(L)-V m (L) satisfies (2.8). To do this, first note 
that from the condition [X(L),L] = 0, L e A, we have [dX(L) ■ L\ L] + [X(L),L'] = 
0, L,L' G A. Therefore, [Y(L),L] = [dX (L) ■ V m (L) , L] = —[X(L),V m (L)] = 
-(m + l)L m [X(L),L] = 0, for all L e A. On the other hand, it follows from 
dX(L) ■ [V, L] = [L', X(L)], L, V e A, that 

(*) (d 2 X(L) ■ L')([L", L]) + dX(L) ■ [L", L'\ = [L" , dX(L) ■ L% L, L', L" e A. 

Consequently, for all L, V e A, we have 

dY(L) ■ [L', L] = (d 2 X(L) ■ [L\ L]){V m {L)) + dX(L) • ((m + l)L m [L\ L}) 

= (d 2 X(L) ■ V m (L))([L', L]) + dX(L) ■ [L', V m (L)} 

= [L',Y(L)], by (*). 

□ 

Remark 4-6. For the vector fields Z n in (2.9), we have in particular the relations 
L Vm Z n = nZ m+n , m > -1, n > 1. 

If we now combine Theorem 4.5 and Theorem 4.2, the nature of the vector fields 
V m is now revealed. 

Corollary 4.7. L\V n = 0, n > 0, L^V- X = 0. □ 

5. Virasoro Action on the Bivector Fields and Compatibility of the 
Structures. 

The goal of this section is to prove Theorem 3.2. To do this, we consider the 
action of the vector fields V m on the bivector fields 7r n corresponding to {■, -}( n ), 
n > -1. 
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Theorem 5.1. Ly m ix n = (n — ra)Tx m+n , m,n> — 1. 

As indicated in Section 3, this result is the key in proving Theorem 3.2. The 
demonstration of Theorem 5.1 is quite tedious, so we break it up into several steps. 
First, note that from the property of the Lie derivative, we have 

(5.2) L Vm 7r n (dF, dH) = V m {F, H} (n) - {V m F, H} (n) - {F, V m H} (n) . 

Using the expressions for {•, and V m , we obtain the identities in the next two 
lemmas. We shall omit the rather lengthy computations. 

Lemma 5.3. 

V m {F,H} {n) (L) = (V m (L), [R(L n+1 dF(L)), dH(L)}) 

+ (L, [R(L n+1 dF(L)), d 2 H(L) ■ V m (L)}) + (n+ 1)(L, [R(L m+n+1 dF(L)), dH{L)]) 

+ (L, [R(L n+1 d 2 F (L) ■ V m (L)), dH(L)}) - (F <- H) 

where (F <-> H) denote terms obtained from previous ones by switching F and H . 

□ 

Lemma 5.4. 

{V m F, H} {n) (L) + {F, V m H} {n) (L) = (L, [R{L n+1 d 2 F{L) ■ V m (L)) , 
dH{L)] + [d 2 F{L) ■ V m {L),R{L n+1 dH{L))]) + (m + 1)(L, [R(L m+n+1 dF '(L)) , dH(L)] 
+ [L m dF{L), R(L n+1 dH(L))]) - (F H) . 

□ 

Proof of Theorem 5.1. By combining the expressions in Lemma 5.3 and Lemma 5.4 
according to (5.2), it is clear that terms involving second derivatives cancel out, and 
we obtain 

(*) L Vm * n (L)(X u X 2 ) 

= (V m (L), [R(L n+1 X 1 ),X 2 ]) + (n - m)(L, [R{L m+n+1 X x ) , X 2 ]) 
- (m + 1)(L, [L m X u R(L n+1 X 2 )]) -(1^2) 
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where X\ = dF(L),X 2 = dH(L). Now, by repeated application of the derivation 
property, the commutativity of multiplication and its symmetry with respect to 
(■, •), we have 

(V m (L),[R(L n+1 X 1 ),X 2 }) - (1 <- 2) 

= (L, [R(L n+1 X 1 ), L m X 2 }) - (LX 2 , [R(L n+1 Xi), L m \) - (1 <- 2) 
= (L, [R(L n+1 Xi), L m X 2 }) - m(L m X 2 , [R(L n+1 X 1 ), L}) - (1 <- 2) 
= (m + 1)(L, [i2(L n+1 X!), L m X 2 ]) - (1 <- 2) 
= (m + R(L n+1 X 2 )\) - (1 <- 2) . 

If we substitute this in (*), the result follows. □ 

Remark 5.5. In the case of noncommutative, associative algebra, relations similar 
to the ones in Theorem 5.1 were obtained in [LP2] for the three structures there. 

Corollary 5.6. [7r m ,7r n ] s = ^[V n+1 , [?r_i, TT m } s }s + m ~+ 2 ~ 1 ^m+n+i]s for 
m,n> —1. 

Proof. From Theorem 5.1 and the graded Jacobi identity for the Schouten bracket, 
it follows that 



1 



Km, [V n +1, 7T-l]s]s 



n + 2 

= — 7J [K+l, k-l, 7Tm]s]s H — k-l, [K+l, 7Tm]s]s 

[ jyx — JX — 1 
= — K Wn+U k-l, 7Tm]s]s H ~ ~ [tT-1, ^m+n+l\S 

n + 2 n+2 

□ 

Remark 5. 7. The formulation of Corollary 5.6 is motivated by similar considerations 
in [ AvM] . 

Proof of Theorem 3.2. If we set m = — 1 in the identity in Corollary 5.5, we find 
[7T_i,7r n ] s = - 2 ( n 1 + 2) i v n+i, [t-i,t-i]s]s = 0, V 7i > -1, as tt_i is the bivector 
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field for the Lie-Poisson structure {•, From the same identity, it now follows 

that [ir m ,iT n ]s = 0, V m, n > —1. Hence the brackets {•, -}( n ) define compatible 
Poisson structures on A. Finally, the assertion in part (c) follows from Theorem 
5.1. □ 



6. Some Examples. 

In this section, we look at some concrete examples of partial differential equations 
which can be realized as Lax equations on Poisson algebras. In each case, we de- 
scribe the multi-Hamiltonian formalism which follows from our universal construc- 
tion in Sect. 3. The reader should note that in these applications, we are dealing 
with Lax operators which form submanifolds of the full Poisson subalgebras under 
consideration. Although these submanifolds of Lax operators are invariant under 
the dynamics of the associated Lax equations, however, they are not automatically 
Poisson submanifolds of the brackets which arise from the general construction in 
Sect. 3. For this reason, there are two kinds of situations in the examples which 
follows. In the happy case where the sub manifold M. of Lax operators does form a 
Poisson submanifold of (A, {•, -}( n )), there is of course an induced structure on M 
which can be obtained by simple restriction of {•, -}( n ) to M. On the other hand, 
when M. is not a Poisson submanifold of (A, {•, -}( n )), the reader will see that the 
geometry in each case warrants the application of Dirac reduction, i.e. reduction 
with constraints [D, KO, MR]. Thus in this latter case, the brackets which arise 
from the construction in Sect. 3 serve as the starting point of a reduction process 
from which the constrained brackets on M. are computed. 

In the following, we shall rescale the expression for {•, -}( n ) by the factor tj. 

1° The Benny hierachy. 

The Benny equations in nonlinear waves [B] (we shall consider the simplest case 
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here) are given by the quasi-linear system 

U \ ( Uq l \ ( Uq 



U-i J t \ U-i U J \ U-i 



(6.1) 

We shall deal with the case where u ,U-i are smooth functions on the circle S 1 



Following [G-KR], introduce the algebra A of Laurent polynomials in A, having 
the form 

(6.2) u(x, A) = ^ u i (x)X i , 

i 

where the coefficients Ui are smooth functions on the circle S 1 . With the well-known 
Lie-bracket defined by 

du dv du dv 

(6 - 3) [uM - 1 = dXd- x -d- x dx^ u ' veA ' 

it is clear that (A, [•, -]_x) is a Poisson algebra. In [G-KR], the Benny equations 
are rewritten as a Lax equation in this Poisson algebra. Indeed, (6.1) is equivalent 
to 



, \ dL 
<6 ' 4 » M 



n I ' /.* I l 



-1 

where the Lax operator L is an element of the Benny manifold 

(6.5) A^Benny = {L E A \ L(x, A) = A + U (x) + U- 1 (x)X~ 1 } 

and the r-matrix R is the one associated with the direct sum decomposition 

(6.6) A = A^@A^ Q 
into subalgebras 

(6.7a) A-^x = < u E A | u(x, A) = ^ u i (x)X i 
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(6.7b) A^o = I u G A | u(x, A) = ^ Ui(x)\ l 

In view of the representation in (6.4), the quasi-linear system (6.1) is only a member 
of a hierachy of Lax equations on M. Benny, and this is what we call the Benny hier- 
achy. Note that the Poisson algebra introduced above admits the trace functional 

(6.8) tr_i-u = J U-i(x)dx , uEA 

(here and below we integrate over S 1 ) 
which satisfies the important property 

(6.9) tr_ 1 [u,v]=Q, u,veA. 

Therefore, we can equip A with a non-degenerate ad-invariant pairing (•, •)_].: 

(6.10) (u, v)-i = tr_i(m;) , u,v<eA. 

Thus we have all the ingredients which are required for the application of Theorem 
3.2. Consequently, we have a family of Poisson structures {•, n > — 1 5 on A. 

It is easy to check that M. Benny is a Poisson submanifold of (A, {■, -}(_i)). There- 
fore, the induced structure on M. Benny provides the first Poisson structure for the 
equations in the Benny hierachy [G-KR]. Using u = (iio,U-i) as coordinates on 
■M Benny, the associated Hamiltonian operator is given explicitly by 

(6.11) B ( . I)( «) = (» D Q ), 

which is apparently well-known to people working in other frameworks (see, for 
example, [DN]). Clearly, this first structure is degenerate, with Casimirs given by 
Ci(u) = f uo(x)dx and C*2(w) = / U-i(x)dx. 

Remark 6.12. One of the advantages in formulating the Benny equations as a Lax 
equation on A is that it automatically suggests a method of solution, namely, via a 
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factorization problem on a symplectic diffeomorphism group. The analytic details, 
however, are nontrivial. 

We now turn to the higher structures. Here, it is easy to see that .MBenny is not 
a Poisson submanifold of any of the brackets {•,•}(«), n ^ 0. However, we shall see 
that we can apply Dirac reduction to {•, -}( n ) with appropriate constraints to obtain 
the higher structures on Alsenny We shall illustrate the procedure for n = and 
n = 1, thereby obtaining the second and third Poisson structures on .MBenny 

For n = 0, the Hamiltonian vector field generated by H is of the form 

(6.13) xP(u)=uTl^ 2 ([dH(u),u]- 1 ) - [n^M# (w)),^ 

= [U^udHiu)), u]_ t - ull>_i([dfl-(u), . 

If L G .MBenny, it follows from this formula that the highest order term of X^\l) 
in A is A , while the lowest order is in A -2 . Using u = (uq, it_i, u_ 2 ) as coordinates 
on the submanifold {X + uo(x) + u-i(x)X~ 1 + tt_ 2 (x)A~ 2 G A}, the operator which 
gives X^(L) can be computed explicitly: 

/ D u D + u 0x U-iD + u-i x 

(6.14) u D 2u-!D + u- lx 
\u-iD — u 2 i 1 D — u-\U-\ 2 

Therefore, we can apply Dirac reduction with constraint -u_ 2 = to obtain the 
second structure on .MBenny: 

(6.15) 

2 u \ n ^f0 l\ fO 0\ 

u 2u_ 1 ) D+ {0 o) U °* + {l 0) U -^ 

Note that this second structure is of hydrodynamic type [DN] because the associated 
Hamiltonian operator is of the form 

(6.16) Bl\u)=^(u)D + ^(u)u k x . 
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In this case, the metric which defines the structure (6.15) is non-degenerate where 

(6.17) A = ug-4u_i^0. 

For n = 1, i.e. for the bracket {•, -}(i), we have a similar formula for the Hamil- 
tonian vector field 

(6.18) x£\u) = u 2 U^_ 2 ([dH(u) 1 u]_ 1 ) - [U^dHiu)),^^ 

= [H^ 1 {u 2 dH{u)) 1 u] _ 1 - ^U^ddHiu), . 

This time, the highest order term of X^\l) (L e -MBenny) in A is still A , but the 
lowest order term is in A -3 . Therefore, in the coordinates u = (uo, it_i, U-2, W-3), 
the operator which gives Xjp(L) is given by 
(6.19) 

2u D+u 0x (uq+3u-i)D+2uoUox+2u-1x 2uqU-iD+2uoU-i x +2u-iUo x u 2 _ 1 D+2u_ 1 u_ 1x 

(uq+3u-!)D+U- 1x 4u u^ 1 D+12u^ 1 u 0x +2u u- 1x u 2 _ 1 D+2u- 1 U- lx 

2 2 2 3 2 

2uqU-iD u_ ± D —2u( S u_ 1 D — 2uqU-\U-\ x —u_ 1 uq x —u_ 1 D—2u_ 1 U—i i 

u 2 _ x D -u 3 _ 1 D-u 2 _ 1 U- lx 

To obtain the structure on -Meenny, we have to use Dirac reduction with the con- 
straints U-2 = 0, us = 0. Accordingly, we have to invert the lower 2x2 block of 

(6.19) : 

/ — 2uqu 2 _ 1 D — 2uqU-\U-\ x — u^uqx — u^_ x D — 2u 2 _ 1 U-\ 
\ -u 3 _ l D - u^iU-ix 

(6.20) / 'D- 1 1 



u~ 



1 



1 JJ — 1 1 1 JJ — 1 y up — 1 .} 

u^T u_i wTi" v - 2 Li 



-l "-i 



Hence the Hamiltonian operator of the third structure is given by 
(6.21) 

fy\ — ( 2u D+u 0x (uq+3u-i)D+2u u 0x +2u-i x 



2ti li-l-D + 2li_llio a ;+2lloW-la: M_ 1 -D + 2M_lM_l a! \ / u _l "_1 

iii 1 D+2i i _ 1 ii_ la: y I -^-D~ x -^ g— D _1 -^2 



( 



—i—D 



-l i 



/ 2u ii_iD u_ 1 D \ 
u 2 _ 1 D J 



4uq u +4u 
Kq+4m_i 4upu_ 



D+( 22 U Z) u 0x +(° 22 l o ) u_ la! 
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Again, this corresponds to a bracket of hydrodynamic type and the non- degeneracy 
of the metric is characterized by the same condition in (6.17). 

Remark 6.22. (a) Alternatively, on the symplectic leaves of the first structure 
defined by the conditions J uo(x)dx = const, J U-\{x)dx = const, B-i is invertible 
and therefore we can compute the recursion operator 

(6.23) n = B Q BZ 1 1 = ( 9 U0+ 4 U ^ D " 2 ) 

\2lL-! + U- lx D u J 

^From this, we can check that B\ = TZBq. 

(b) In principle, one can compute all higher structures explicitly by applying 
Dirac reduction to {•, -}( n ) or by using the recursion operator 1Z, but the calculations 
are quite formidable and we do not know if there exists an efficient way to do this. 

2° The dispersionless Toda lattice hierachy. 

Let A be the algebra introduced in Example 1°, but now we equip it with the 
following Lie bracket 

(6 ' 24 > ^- x {axlT x -^Tx) • "' V€A - 

Then (A, [•, -] ) is also a Poisson algebra. The dispersionless Toda lattice hierachy 
is defined by the Lax equations 

(6.25) ^ = [U i (L n ),L] = -[U l (L n ),L} , n = l,2,... 
where the Lax operator L is an element of the manifold 

(6.26) -MdToda = {L G A | L(x, A) = Ul (x)X + u (x) + ui(x)X~ 1 } 

and LTe, LT[ are the projection operators relative to the direct sum decomposition 



(6.27) 



A = t®i 
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into subalgebras 

(6.28a) t = 



{u e A | u(x, A) = ^2 Ui(x)(X l - A l ) 
i>0 



(6.28a) (= I u e A \ u(x, A) = ^u^x)X 

[ i<0 

When n = 1, the corresponding Lax equation 

jt ( u ot = 4mu lx 

(6.29) - F = [ML),L] ^\ 

{ U lt = U!U 0x 

These are the dispersionless Toda lattice equations and can be obtained from the 
periodic Toda lattice ODE system 

(6-30) ^ = 2(bl-bl_ 1 ), ^=b k (a k+1 -a k ) 

by taking a continuum (or long wave) limit. 

The Poisson algebra (A, [•, -] ) also has all the ingredients needed for the con- 
struction in Theorem 3.2. In this case, the invariant trace is of the form 

(6.31) trou = J uo(x)dx , uEA 

which gives rise to the non-degenerate ad-invariant pairing (•, -)q: 

(6.32) (u, v)o = tro(m>) , u,v<EA. 
As the r-matrix for the equations in (6.25) is given by 

(6.33) R = U t -U { , 

it follows from (3.3) that the Hamiltonian vector field generated by H in the struc- 
ture {•, -}( n ) is °f the form 

(6.34) X^\u) = [n t {u n+1 dH(u)),u] Q - u n+1 Ii\{[dH(u),u} ) 

= u n+1 IlU{dH(u),u]o) - [U l {u n+1 dH(u)),u] . 
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Using this formula, we can now check that MdToda is a Poisson submanifold of 
(A, {•, -}( n )) on ly f° r n = — 1, 0, and 1. Accordingly, the induced structures on 
A^dToda provide the first, second and third Poisson structures for the equations 
in the dispersionless Toda lattice hierachy. Using u = (uq,ui) as coordinates on 
-MdToda, the Hamiltonian operator of the first structure is given explicitly by 

^ *_,(„) = (» ;) B+ (° »-). 

Clearly, the associated Hamiltonian vector fields preserve the sign of u\. Therefore, 
B-i(u) restricts to a structure on 

(6.36) M^ Toda = {L e A | L(x, A) = u 1 (x)X + u (x) + wi(x)A _1 , m{x) > 0} 

whose symplectic leaves are the level sets of the Casimirs / uo(x)dx, J \nui(x)dx. 
Finally, we note that S_i(w) is obviously of hydrodynamic type and the corrre- 
sponding metric is non-degenerate on M-^ Toda as det (g 1 -*) = —u\. 

Remark 3.37. Note that the equations in the dispersionless Toda lattice hierachy 
are (strictly) hyperbolic in A^^" Toda and we can take 

(6.38) w\{u) = uq — 2u\ , W2(u) = u$ + 2u\ 

as the Riemann invariants. We shall not give the proof as the reader can easily 
supply the details. 

As for the second and third Poisson structures on MdToda.-, direct calculation 
shows the corresponding Hamiltonian operators have the form 

(6.40) 

/ 8u uj 4itf + ulux \ ( 4uf \ / 8w wi ul + 8uf \ 

V Au\ + u 2 oUl 2u q u\ J + { 2u oUl u\ J UOx+ y Au\ 2u QUl J Ulx ' 
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These structures also restrict to M. J Toda , and are obviously of hydrodynamic type. 
But in contrast to the first structure, the metrics associated with Bq(u) and B\{u) 
are non-degenerate only on a subset of M. ^ Toda , characterized by the condition 

(6.41) w 1 {u)w 2 {u) 0, 

where w\ (u) , w 2 (u) are the Riemann invariants in (6.38). 

In order to compute the higher structures, we have to invoke Dirac reduction, as 
in the last example. Here, we shall do this for the fourth structure as it presents 
new features which are also shared by all higher structures. First of all, we check 
that for L e ■M'^ rod& , we have X^(L) e /mil*, and the highest order term in A 
is A 2 . Then we write down the operator which gives X^ (L) using the coordinates 
u = (uo,ui,u 2 ) on the submanifold where X H '{L) lies: 
(6.42) 

(12u\+12u^u\)D+&{u\+v%u\) x {u^u\+12uou\)D+Qu\uo x +2Auou\ui x +u^ > ui x 2u\D+%u\u\ x 
(uqUi+12uou\)D+u\ (uQ+6u uf) x (Au\+3u^u\)D-\-3u u\u 0x -\-(%u\-\-3u'lu 1 )u lx u\uq x 

2u\D -UjUoi — u\D— 2u\ui a 

Finally, we invoke Dirac reduction with constraint u 2 = to compute the Hamil- 

tonian operator of the fourth structure, and the result is 

(6.43) 

r ( \ - ( 16u i + 12 «o«i 12uo«i + ulii! \ ( 12u u\ Au\ \ 
B2 [U) ~ { 12u ul + u 3 oUl 4uj + 3u 2 ul J + \ Su 2 oUl + 8uf Zu Q u\ ) U ° x 

32uf + 12ulm 24u u 2 + itg A 
12u uj 8uf + 3u 2 oUl J Ulx 



lQu\U\ x D 1 u\Ui x 4u\U\ x D u^uq 

4u X U^ x D~ X UxUx x UiU 0x D~ 1 UiU a 



Thus, B 2 (u) has a nonlocal tail, and provides an example of a class of nonlocal 
Hamiltonian operators of the form 

N 

(6.44) 5* («) = g ij D + 6?«* + ^OOW^^OOM • 

a=l 

In the case where det((7 lJ ) ^ 0, we note that the geometric root of such structures 
was discussed in [F] and applied to the chromatography equations. At this point, 
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the reader can check that the subset of M^ Toda where the metric associated with 
B2(u) is on-degenerate is likewise defined by (6.41). Also, on the symplectic leaves of 
the first structure where f uo(x)dx = const and f \nu\{x)dx = const, the recursion 
operator R = BqBz\ exists and it is not hard to show that Bi = TZbo and B 2 = 
U 2 B . 

Remark 6.45. In [DM], the authors considered the dispersionless Toda lattice equa- 
tions with boundary conditions tti(O) = 0, Wi(l) = 0. We remark that the multi- 
Hamiltonian formalism of this problem can also be obtained in a similar fashion. 
Indeed, the only major change one has to make here is to replace the algebra above 
by the algebra of Laurent polynomials in A, having the form u(x, A) = J2 i Ui(x)X z , 
where the coefficients itj are smooth functions on i" = [0, 1] satisfying the additional 
conditions %(0) = itj(l) = 0, j ^ 0. Otherwise, everything goes through just the 
same as before. In particular, the formula for the Hamiltonian operators of the first 
four structures are still those given in (6.35), (6.39), (6.40) and (6.43). 

In the next two examples, we shall consider equations with infinitely many field 
variables. For simplicity of exposition, we shall not get into reduction calculations 
here, only remark that the number of constraints is still finite in each case. 

3° The dispersionless KP hierachy. 

Let A be the algebra of formal Laurent series in A, having the form 

N(u) 

(6.46) u(x,X) = u i( x ) X ^ 

i= — oo 

where the coefficients u\ are smooth functions on 5" 1 = R/Z. Define 

(6.47) [«,»]_, = „, v eA, 

then (A[- 7 -]_i) is a Poisson algebra. The (extended) dispersionless KP (dKP) hier- 
achy is defined by the equations 

(6.48) ^ = [n> (L n ), L]^ = -[U<_ 1 (L n ),L}_ 1 , n = 1,2,..., 
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where the Lax operator is an element of the (extended) dKP manifold 

(6.49) M dK p =<LeA \ L(x, A) = A + u * 0*0*4 > 

^ i= — oo J 

and LT>o,n<_i are projection operators relative to the decomposition 

(6.50) A = A> ®A<- 1 
into subalgebras 

(6.51a) A> = <j u E A \ u(x, A) = ^ it* (a;) A* 

[ i>0 

(6.51b) A<_i = \u e A | u(x, A) = Y u i( x ) xi 

(. i= — oo 

In the standard form of the dKP equations [TT], the coefficient u$ = 0, but we shall 
not get into reduction calculations here. 

For the Poisson algebra (A, [•, the invariant trace is defined by 

(6.52) tr_iit = J u-i(x)dx , u£A, 
and we have the non-degenerate ad-invariant pairing (•, -)_i: 

(6.53) (u, = ti-i(uv) , u,v G A . 
So again we can invoke Theorem 3. 2, using the r-matrix 

(6.54) i? = n> -n<_! 

in this case to obtain the corresponding brackets {•, -}( n )) n ^ — 1- Here, it is easy 
to check that MdkP is a Poisson submanifold of (A, {•, -}( n )) only for n = —1,0. 
Therefore, the induced structures on MdKP provide the first and second Hamil- 
tonian structures for the equations in the hierachy. For the bracket {•, -}(i), the 
slightly larger manifold jtt e A | u(x, A) = J2l=-oo Ui(x)X 1 ^ is a Poisson submani- 
fold. Hence the third structure on M.dKP can be computed using Dirac reduction 
with constraint u\ = 1. We shall leave the details to the interested reader. 
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4° The dispersionless modified KP and the dispersionless Dym hierachy. 

Let (A, [•, -]_i) be the Poisson algebra in Example 3°, with the same invariant 
pairing (•, Consider the decomposition 

(6.55) A = A>k © A<k-i , k > 

with associated projection operators Tl>k and LT<fc_i, where 

(6.56a) A> k = < u e A \ u(x, A) = ^ Ui(x)X l 

[ i>fc 

{fc-i 
«gA|«(x,A)= ^ u i (x)X i 
i= — oo 

Clearly, A> k is a subalgebra of (A, [•, ■]_{) for all fc. On the other hand, simple 
verification shows that A<k-i is a subalgebra of (A, [-,-]-i) only for k = 0,1,2. 
Therefore, among the direct sum decompositions in (6.55), only the three cases 
k = 0,1, and 2 lead to r-matrices, and the case k = has already appeared in 
Example 3°. We now consider the other two cases, with Lax equations 
(6.57) 

^ = \n> k (L n ), L]_ t = - [n< fc _!(L"), L]_ x , n = 1, 2, . . . , ; k = 1, 2, . 

For k = 1 and L e MdKP, the equations in (6.57) constitute the dispersionless 
modified KP hierachy. For k = 2, we obtain the dispersionless Dym hierachy when 
the Lax operator L is from the submanifold 

(6.58) M dDym = \l e A | L(x, A) = ^ ^(x)A l | . 

These hierachies are the semi-classical limit of the modified KP and the Dym hi- 
erachies in [ANPV,KO]. For the dmKP hierachy, with r-matrix given by R = 
n>i — n< , the manifold of Lax operators is a Poisson submanifold of the asso- 
ciated brackets {•, -}( n ) f° r k = —1,0,1. Hence the induced structures on AidKP 
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provide the first three Poisson structures for the Hamiltonian description of dmKP. 
The higher structures, on the other hand, have to be computed using Dirac re- 
duction. For the dispersionless Dym hierachy, the situation is even better, for in 
this case the first five Poisson structures on MdDym are obtained from the brackets 
{'i '}(n) ( — 1 < ^ < 3) associated with R = n> 2 — II<i by simple restriction. Again, 
the passage from {•, -}( n ) (n > 4) to the higher structures require the application of 
Dirac reduction. 
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